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We study the ground state fidelity, fidelity susceptibility and quench dynamics of the extended 
quantum compass model in a transverse field. This model reveals a rich phase diagram which 
includes several critical surfaces depending on exchange couplings. We present a characterization of 
quantum phase transitions in terms of the ground state fidelity between two ground states obtained 
for two different values of external parameters. However, we derive scaling relations describing the 
singular behavior of fidelity susceptibility in the quantum critical surfaces. Moreover, we study the 
time evolution of the system after a critical quantum quench using the Loschmidt. We find that the 
revival times of Loschmidt echo are given by Trev = Nl2vmax, where N is the size of the system 
and Umax is the maximum of lower bound group velocity of quasi-particles. Although the fidelity 
susceptibility shows the same exponent in all critical surfaces, the structure of the revivals after 
critical quantum quenches displays two different regimes reflecting different equilibration dynamics. 

PACS numbers: 03.65.Yz, 05.30.-d, 75.10.Pq 


I. INTRODUCTION 

In the last few years a big effort has been assigned to the analysis of quantum phase transitions (QPTs) from 
the perspective of quantum information [IHE]- Entanglement and fidelity have been accepted as new notions to 
characterize quantum phase transitions. Entanglement, referring to quantum correlations between subsystems is a 
good indicator of quantum phase transitions, because the correlation length diverges at the quantum critical points 
m- The fidelity which is a measure of distance between quantum states, could be also a nice tool to study the 
drastic change in the ground states in quantum phase transitions [13]. 

On the other hand, recent advances in the studies of ultra-cold atoms trapped in optical lattices introduced a 
new tool to simulate the dynamics of interacting quantum many-body systems in non-equilibrium strongly correlated 
quantum systems [Til ITS] . These new opportunities are concentrated by the pioneering experiments on tunable 
Mott insulator to superfluid quantum phase transitions, observed by utilization of the optical lattice potential in 
three-dimensional 3D El and ID [TB] systems. 

Recently the investigation of non-equilibrium properties of closed quantum systems have been getting a lot of 
attention for several reasons. Specifically it has been applied to quantum information in which decoherence and 
entanglement dynamics play an essential role. Erom the theoretical point of view, it is very important to understand 
the notion of universality for a system away from equilibrium, where the traditional concepts of phase, fixed point and 
renormalization group fail. Driving a system out of equilibrium is done in many ways. Most of the attention has been 
focused on quantum quenches na, namely, sudden changes of the external parameters of the Hamiltonian controlling 
the unitary evolution of the closed system. One of the conventional methods of understanding the dynamics of a system 
after a quench is the Loschmidt echo (LE), which is a benchmark of the partial or full reappearance of the original 
state as a function of time [HESI. The LE is defined as follows: if a quantum state \ip) evolves with two Hamiltonians 
H and H', respectively, the LE is the measure of the overlap given by LE = |('0| exp(iiLT) exp(—If the 
system admits the ground state, LE is a dynamical version of the ground state fidelity. Recently, the time behavior 
of the LE has been studied in well-known models, in particular the XY spin chain [iniisi] and cluster XY chain 
[22| . High values of the LE mean that the system is approaching the initial state. Typically, the LE will decay 
exponentially at first and then start oscillating around an average value |20j . If the system is finite, the time evolution 
is quasi-periodic, forcing the system arbitrarily close to the initial state for long enough times. The system will show 
revivals, i.e., times when the value of the LE is greater than the average value. The structure of these revivals may 
be greatly influenced by criticality [20] . In this work, the phase diagram and universality of the one-dimensional 
extended quantum compass model (EQCM) [23H25] in a transverse filed [26| |27] will be studied by means of the 
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fidelity of the ground state and fidelity susceptibility. This inhomogeneous model covers a group of well-known spin 
models as its special cases and shows a rich phase diagram. It should be mentioned that the phase diagram [25], 
quantum correlation [28] . bipartite entanglement |29j and fidelity |30j . of this model have been studied numerically 
using the exact diagonalization method and infinite time-evolving block decimation |32] . This study is an important 
addition to the literature as to the best of our knowledge, the quenches of a Hamiltonian which also undergoes a very 
unique type of quantum phase transition [23l|26][27], has not been investigated before. However, it will be instructive 
to study the quenches near quantum critical points because of the expected universality of the response of the system, 
and thus the possibility of using the quench dynamics as a nonequilibrium probe of phase transitions. Then, we will 
explore the quenching dynamics of this model, when the transverse field or the exchange couplings is quenched. 


II. THE HAMILTONIAN AND ITS EXACT SOLUTION 

Consider the Hamiltonian 


N’ 

n—1 


^2n) 


( 1 ) 


where Ji and J 2 are the odd bonds exchange couplings, Li and L 2 are the even bond exchange couplings and N = 2N' 
is the number of spins. It should be pointed out that although periodic and antiperiodic boundary conditions differ 
in 0(1/N) terms, this difference usually does not affect the phase diagram or other quantities in the thermodynamic 
limit. However, although it can be important in the LE that is typically exponentially small in N and the boundary 
conditions could have a dramatic effect in the case of the critical quench |3Ij . but for simplicity we assume periodic 
boundary conditions. This model embraces a group of the other familiar spin models as its special cases, such as 
the quantum Ising model in a transverse field for J 2 = L 2 = 0, the transverse field XY model for Ji = Li and 
J 2 = L 2 , and the transverse field XX model for Ji = J 2 = Li = L 2 . The above Hamiltonian [Eq. 0 ] can be exactly 
diagonalized by standard Jordan-Wigner transformation [571133] as dehned below. 
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which transforms spins into fermion operators Cj. 

The crucial step is to define independent Majorana fermions [331|3S] at site ri, = C 2 „_i and = C 2 n- This can 
be regarded as quasiparticles’ spin or as splitting the chain into bi-atomic elementary cells [51] . 

Substituting for crj, trj and cr| (j = 2n, 2n— I) in terms of Majorana fermions followed by a Fourier transformation, 
Hamiltonian Eq. 0 (apart from an additive constant), can be written as 


H+ = 


^ [jcl^cP_\ + Lcl^4 + 2h(cf c« + + h.c. 


where J — {Ji — J 2 ) — (-^1 ~ 5 ^ — (<^1 J 2 ) (-^1 and k — (j — 3, * * * , N' — 1). 

By grouping together terms with k and —k, the Hamiltonian is transformed into a sum of independent terms acting 
in the 4-dimensional Hilbert spaces generated by k and —k {H'^ — 0fc>o other word [Hk, Hk'] = 0 in which 

= (JcfVI, + Lcf4 - J’cy_, - L‘44< + Jc%4} + - J‘c« (2) 

+ 2Mcf 4 + ey, + +iV-P- 

Hamiltonian Eq. 0 can be written in the diagonal block form 
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where rj) = (c«, c!.\, c^, and 
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FIG. 1: (Color online) Three-dimension of the ground state fidelity of quantum compass model in a transverse field as a function 
of magnetic field h and Hamiltonian parameter J 2 . The parameters set as (a) Ji = Li = 1, L 2 = 0, N = 100 and Sh — 0.001, 
(b) Ji = —1, Li = 1 L 2 = 0, = 100 and Sh = 0.001. 


By using the element of the new vector FJ! = ( 7 ^, 7 ®^,, 7 ^, 7 !!.^,) which could be described by unitary transformation 
r'j, = UkFk (see Appendix A), the matrix A(k) can be diagonalized easily and we find the Hamiltonian Eq. (|^ in a 
diagonal form. 


h = E [e|( 7|'72 - i) + ££(7*71 - 5 ) 

k 

where = \/a -I- y/b and = yja — y/b, in which 


(4) 


a = ih^ + I Jp -b |Lp, b = (16/1^ -b 2| J^lLp -b 

The ground state {Eq) and the first excited state (Ei) energies are obtained from Eq. (§, 

k k 

It is straightforward to show that the energy gap vanishes at ho = yj (Ji -b L 2 ){J 2 + Ti) and = 
yj {Ji + T 2 )( J 2 — Li) in the thermodynamic limit. 

So, the quantum phase transition which could be driven by the transverse-field depending on exchange couplings, 
occurs at ho and hj^- 

By a rather lengthy calculation on the unitary transformation we can obtain the whole spectrum and the eigenstate 
of the Hamiltonian which have been written in the vacuum fcth mode of c| and c^, 


IV’m) 


= nKI0) + < + < cfc^t,|0) + < c'^^A^\0) + vT + cfcflO) 

k 

+vT c-J,cti\ 0 )+v^ 4’=lU’<ibo>l. 


(5) 


where \tpm) {m = 0, • • • ,7) is the eigenstate of the Hamiltonian with corresponding eigenvalue and vj, {j = 
1, • • • , 8 ) is functions of the coupling constants (see Appendix B). 
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FIG. 2: (Color online.) Contonr plot of the fidelity of quantum compass model in a transverse field. The parameters set 
as Li = 1, L 2 = 0, = 400 and (a) h = 0.1, (b) h = 0.5. The light blue line represents ho = \/J 2 + 1 critical line and 

Att = y/J 2 — 1 critical line has been shown by the red lines. 


III. FIDELITY AND PHASE DIAGRAM 

The phase diagram can be examined by considering the fidelity and fidelity susceptibility introduced in Ref. m- 
As mentioned in introduction, the fidelity of ground state is defined by the overlap between the two ground state 
wave functions at different parameter values and is defined as 

F{h,Sh) = \{Mh)\Mh + Sh))\, (6) 

where |'0o(^)) is a ground state wave function of the many-body Hamiltonian describing the system exposed to an 
external magnetic field h, while Sh is a small deviation from h. The main idea is that near a QPT point there is a 
sharp enhancement in the degree of distinguishability between two ground states, corresponding to different values of 
the parameter space which defines the Hamiltonian. This distinguishability can be determined by the fidelity, which 
for pure states simplify to the amplitude of inner product or overlap. The behavior can be ascribed to a sudden 
change in the structure of the ground state of the system across the quantum phase transition. Therefore, one expects 
that fidelity has a minimum at the critical point and it should contain all the information that describes QPTs and 
topological order. The drop of fidelity specify not only the position of the critical point, but also universal information 
about the transition given by the critical exponent u which is the correlation length exponent associated with the 
QCPs [13]. The response of the hdelity after an infinitesimal change of the external parameter up to second order 
reads 


F{h, Sh) = 1 — 


Sh^ 

-^XF, 


where the fidelity susceptibility {xf) is defined by [36ll37] 


XF = {dhtjJo{h)\dhtpo{h)) - {dh'(po{h)\tpo{h)){tpo{h)\dhtjJo{h)). (7) 

If there exist more than one external parameter, this result could be generalized to the so-called quantum geometric 
tensor [I31I3MI. 

The fidelity of this model has been calculated using Eqs. (§ and (§. Three-dimensional panorama of the ground 
sate fidelity of the model has been depicted in Figs. (a), (b) versus the magnetic field and J 2 for Ji = 1 and 
Ji = —1, where we have set N = 100 and Sh = 0.001. Obviously, there is a sudden drop in the ground state fidelity 
at the QPTs lines. 
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FIG. 3: (Color online.) Contour plot of the fidelity of quantum compass model in a transverse field. The parameters set as 
Li = 1, L 2 = 0, = 400 and (a) h = 1, (b) h = 1.5. ho = -\/J 2 + 1 and hr: = — 1 critical lines have been shown by the 

light blue and the red lines respectively. 


TABLE I: The critical exponents 6, t, v. The scaling behavior of the susceptibility in the vicinity of the critical points. 


Exchange couplings 

Critical points 

9 

T 

V 

Ji = 1, J 2 = 0.44 

he — ho —' 1.2 

1.78 

2.02 

1.00 

Ji = 1, J 2 = 2 

0 

t-H 

II 

II 

u 

1.90 

1.97 

1.00 

Ji = 1, J 2 = 2 

11 

0 

II 

0 

1.90 

2.01 

0.985 

Ji = -1, J 2 = 0.51 

0 

II 

II 

u 

2.70 

1.93 

1.01 


In Fig. (a), one observes that the transition lines ho = ■v/J 2 + 1 and = \/— are characterized by two 
assumed lines on the minimum parts of the fidelity surface. However, in Fig. (b) the minimum line on the fidelity 
surface defined by ~ 1 clearly indicates a second order phase transition. By computing the fidelity of the 

one-dimensional extended quantum compass model in a transverse filed, we find the expected critical lines that we 
already discussed in the context of the exact solution and we illustrate the phase diagram of this model in Figs. 
(a), (b),| (a), (b). In Ref. [26] the phase diagram of this model has been investigated by use of the gap analysis and 
universality of derivative of the correlation functions. As previously mentioned, the fidelity and fidelity susceptibility 
of the ground state properties could reflect different zero-temperature regions. This model is always gapful except at 
the critical surfaces where the energy gap disappears. There are four gapped phases in the exchange couplings’ space: 

• Region (I) Ji > 0, 0 < J 2 < J^iJijh): In this region for h < ho the ground state is in the spin-flop phase (Figs. 
0 and ([3|). 

• Region (II) Ji >0, J 2 > J|( Ji, /i): In this case there is antiparallel ordering of spin y component on odd bonds 
((^yy ^ 0)[23l |2ll [26] for h < h^. In this region tuning J 2 forces the system go into a spin-flop phase (region 

(I))- 

• Region (III) Ji < \ Jf{J 2 ,h)\, J 2 > J|(Ji,/i): In this region the ground state is the ferromagnetically polarized 
state (Saturate Ferromagnetic) along the magnetic field (SF). (Figs. (|^ and ^). 

• Region (IV) Ji <0, J 2 < J^iJi^h): In this region the ground state is in the strip antiferromagnetic (SAF) 
phase for h < h^. 
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FIG. 4: (Color online) (a) Fidelity susceptibility per particle of quantum compass model in a transverse field as a functions of 
h for various system size for Ji = 1, J 2 = 0.44. (b) Scaling of the maximum of Xf in terms of system size (A^). Inset: Scaling 
of the position (JiMax) of xf for different-length chains where hMax is the position of maximum in Fig|^(a). 




FIG. 5: (Color online) (a) The finite-size scaling analysis for the case of logarithmic divergence around the maximum point 
(hMax) for Ji = 1, J 2 = 0.44. The fidelity susceptibility, considered as a function of system size and couplings, collapses on a 
single curve for different lattice sizes, (b) Fidelity susceptibility of the model as a function of h for different system sizes for 
Ji = 1, J 2 = 2. 


IV. UNIVERSALITY AND SCALING OF FIDELITY SUSCEPTIBILITY 

It is expected that the fidelity susceptibility (xf) probes the QPTs. Then it will be helpful to study the universality 
and scaling behavior of xf to better understand the properties of the fidelity, and the relation between fidelity and 
quantum criticality. In this section we investigate the scaling behavior of the fidelity susceptibility by the finite size 
scaling approach. In Fig. Q (a) the two dimensional plot of the fidelity susceptibility per particle {xf/N) has been 
shown versus the magnetic field for different system sizes for Ji = 1, J 2 = 0.44. 

Although there is no real divergence for hnite lattice size, but the curves exhibit marked anomalies with height of 
peak increasing with the system size and in thermodynamics limit xf/N diverges as the critical point is touched. 
More information can be obtained when the maximum values of each plot and their positions are examined. As 
it manifests the divergences of xf occurs at he = 1.2 where exactly correspond to the critical point that has been 
obtained using the energy gap analysis {he = ho = ^/lA4). Our investigation manifest the scaling behavior of xf at 
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FIG. 6: (Color online) Three-dimension of the Loschmidt echo as a functions of magnetic field and time, where the Hamiltonian 
parameters set as Ji = 1, J 2 = 1.25, Li = 1, 1/2 = 0 for N = 400 and Sh = 0.0001. 


(a) (b) (c) 



FIG. 7: (Color online.) Variation of ground state LE versus the time t to the critical point Ji = 1, J 2 = 1.25, Li = 1,L2 = 
0,h = 0.5 with N = 400, starting from different values of the coupling constants and magnetic field (a) Ji = 0.8, Ji = 1.2, (b) 
J 2 = 1.2, J 2 = 1.4, (c) h = 0.3 and h = 0.7. 


the maximum point versus N. In this way we have plotted the scaling behavior of XplhMaa: |5(b)) which shows 

a linear behavior of ^"t^{xF\hMai) versus In(iV) with the exponent t ~ 2.02 ± 0.01. 

A more detailed analysis shows that the position of the maximum point (hMax) of xf tends toward the critical 
point like hMax = he + N~^ {6 — 1.78 ± 0.01) which has been plotted in the inset of Fig. |4](b). 

To study the scaling behavior of fidelity susceptibility around the critical points, we perform finite-scaling analysis, 
since the maximum value of xf scales logarithmically. Then, by choosing a proper scaling function and taking into 
account the distance of the maximum of xf from the critical point, it is possible to make all the data for the value of 
xfI^'^I'' as a function of {hMax — he) for different N collapse onto a single curve. The analysis of the finite-size 
scaling is shown in Fig. (a) for several typical lattice sizes. It is clear that the different curves which correspond 
to various system sizes collapse to a single universal curve as expected from the finite size scaling ansatz. Our result 
shows that v = \ exactly corresponds to the correlation length exponent of Ising model in a transverse field {v = 1). 

We have plotted Xf/N for Ji = 1 and J2 = 2 versus h in Fig. (b) for different lattice sizes which shows the 
singular behavior as the size of the system becomes large. As it manifests the divergences of xf occurs at he^ = 1 
and he 2 = ■\/3- The similar analysis shows the scaling behavior of the position of the first and second maximum point 
{hMaxi, hMax 2 ) of XF tends toward the critical point like hMaxi = he^ + [i = 1, 2) and results are presented in 
Table. (I). Moreover, our results show a linear behavior of versus In(A^) with the exponent Ti (see Table I). 

We illustrate the finite-size scaling behaviors of xf around its maximum points. It shows that the fidelity susceptibility 
can be approximately collapsed to a single curve. These results show that all the key ingredients of the finite-size 
scaling are present in these cases too. In these cases scaling is fulfilled with the critical exponent v = 1 (see Table. (I)), 
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FIG. 8: (Color online.) The LE for the quenched quantum compass model to the critical point Ji = 1, J 2 = 1.25, Li = 1, L 2 = 
0, /i = 1.5 with N = 400, starting from different values of parameters (a) Ji = 0.8 and Ji = 1.2, (b) J 2 = 1.1 and J 2 = 1.4, (c) 
h — l.‘i and h — 1.7. 



t t t 


FIG. 9: (Color online.) The LE starting from the SAF and SF phases respectively quenched to the critical point Ji = —1, J 2 = 
0.75, Li = 1,1,2 = 0, h = 0.5 with N = 400, (a) Ji = —1.2 and Ji = —0.8, (b) J 2 = 0.6 and J 2 = 1, (c) h — 0.3 and h — 0.7. 


in agreement with the previous results and the universality hypothesis. The behavior of the fidelity susceptibility has 
been also investigated in other regions. Our calculations show that the non-analytic and scaling behavior of fidelity 
susceptibility are the same as the former results with the same finite size scaling (see Table I). It would be worth to 
mention that the defect density nex{h) can be related to the fidelity susceptibility xf as riex = {^/L'^)hxF{h), where 
L is the linear dimension of a d-dimensional quantum system. The defect density can be defined for any Hamiltonian 
system, for a sudden quench close to a quantum critical point [JD]. In a sudden quench, when a parameter h in 
the Hamiltonian of the system is changed suddenly, the wave function of the system does not have sufficient time to 
evolve. If the system is initially prepared in the ground state for the initial value of the driving parameter, it can 

(a) (b) (c) 




FIG. 10: (Color online.) The LE starting from the different phases to the critical points (a) Ji = 1, J 2 = 1.25,1,1 = 1,L2 = 
0,h = 0.5, (b) Ji = 1, J 2 = 1.25, h — 1.5, (c) Ji = —1, J 2 = 0.75, h = 0.5, for different system sizes. 


























































































9 





FIG. 11: (Color online.) The scaling behavior of the first revival time Trsv for different-length chains for a quenches at 
the different types of critical points (a) Ji = 1, J 2 = 1.25,/i = 0.5, (b) Ji = 1, J 2 = 1.25, Li = 1,1/2 = 0, ft = 1.5, (c) 
Ji = -1, J 2 = 0.75, 0.5 . 


(a) (b) 





FIG. 12: (Color online.) The first derivative of ground state (Eq) and excited state (Ei) with respect to the momentum k 
for quenches at the different critical points (a) Ji = 1, J2 = 1.25, Li = 1,L2 = 0,h — 0.5, (b) Ji = 1, J2 = 1.25, h = 1.5, (c) 
Ji = -1, J 2 = 0.75, h^O.5 0.1, (b) h = 0.5. 


not be in the ground state of the final Hamiltonian. Consequently, there are defects in the final state and its scaling 
follows from the scaling of fidelity susceptibility. 


V. QUENCH DYNAMICS AND LOSCHMIDT ECHO 

Studying the quench dynamics of the systems can be done in several methods. One of the possible scenarios is 
Kibble-Zurek mechanism [41] which estimates the functional dependence of the density of defects on the quenching 
rate in a system crossing a quantum critical point. Moreover, the Landau-Zener formula [42] turns out to be very 
valuable in investigating the transition probability to the excited state of a two-level systems when a parameter in 
the Hamiltonian is changed at a slow and uniform rate. Another formalism which has been applied in this section to 
investigate the time evolution of the ground state after a critical quantum quench is Loschmidt echo. It should be 
mention that, although the Landau-Zener formalism is a powerful tool to study the quench dynamics of systems, it is 
not possible to deal with the Hamiltonian of the model Eq. ([^ as a two-level system. 

As has been mentioned in introduction, a quantum quench is a sudden change in the Hamiltonian of the system. 
The quantum system is originally prepared in the ground state |'0o) of and at time t = 0 the parameters are 

switched to different values The system then evolves unitarily with the quench Hamiltonian according 

to IV'o(O) = where U{t) = exp{—iH{h^‘^^)t). An important quantity describing the time evolution is the 

Loschmidt echo (LE) defined as 

which gives a measure of the distance between the time evolved state \ijjo{t)) and the initial state IV'o)- High 
values of the LE mean that the system is approaching the initial state {LE{t = 0) = 1). The LE typically decays 
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exponentially in a short time T^ei (relaxation time) from 1 to its average value, around which it oscillates. Revivals are 
also visible in the LE as deviations from the average value. The revivals have been defined as a time instances {Trev) 
at which the signal LE{t) differs from the average value. The structure of these revivals could be greatly affected by 
criticality m- Since the model is exactly solvable and the whole spectrum and the eigenstate of the model has been 
obtained, this leads to an exact expression for the LE, and the revival times can be extracted by inspecting its time 
dependence. To obtain the analytical expression for LE, we should express the old ground state in terms of 

the eigenstates of the quench Hamiltonian The time evolution in the EQC model in a transverse field after 

a quantum quench is given by: 


IV’o(O) = - 




7 ^ 


where = \/|(V’m|V'm)l; (m = 0, • • • ,8), and the LE is thus 






E 




( 2 )^ 






i=i 


In the following, the detailed analysis of the LE has been shown for different types of quenches. Three-dimensional 
panorama of the LE has been plotted as a functions of a transverse field and t in Fig. [^for Ji = 1, J 2 = 1.25 and 
system size set as TV = 100. As one can see, LE experiences a sharp decay around the critical points = 0.5 and 
/iQ = 1-5 which exactly correspond to the critical points that have been obtained using the fidelity and gap analysis 
/itt = \/ ( Ji -f L 2 ){J 2 — Li), /iQ = \/ ( Ji -I- L 2 ){J 2 + Li). As it is clear, the structure of revivals are different at the 
critical points and T^eviho)/Tj-evih-n) — 2.5. 

Consider the critical pint Ji = 1, J 2 = 1.25,/i = 0.5. As we see in Fig. [^(b), this point lies between two different 
phases. Increasing (decreasing) Ji and J 2 (magnetic field h) forces the system into the antiparallel ordering of the y 
spin component on even bonds phase and decreasing (increasing) Ji and J 2 (magnetic field h) drive the model to the 
spin-flop phase. In Fig. [^LE has been plotted for a different types of quenches. For this critical point the behavior 
of the LE does not depend on which phase the system is prepared in qualitatively, and T^ev is the same for all the 
quenches. 

In Fig. [^the system is quenched from different ground states corresponding to different phases (saturate ferromag¬ 
netic phase and spin-flop phase) to the critical point Ji = 1, J 2 = 1.25, h = 1.5. The numerical simulations show that 
the details of the evolution are somewhat different, but the structure of the revivals is the same for all the quenches, 
and it oscillates around a relatively high mean value in Fig. (b) in compare with those one in Figs. (a) and|8 (c). 
The behavior of the LE after a quench of Ji and magnetic field (Figs. (a) (c)) are similar to the one obtained for 
critical quenches in the XY model [S]. 

We also quench to the critical point Ji = —1, J 2 = 0.75, h = 0.5 which lies between SAF phase and SF phase (Fig. 
([^). In this case, the details of the LE evolution are the same too, and the only difference is their oscillations around 
a different mean values. These indicate the universality of the revival structure in which the initial state and the 
size of the quench are unimportant. The question which is relevant in the measurement process is how the revivals 
time (Trev) can be derived in the LE? We can address this question for the class of quasifree Fermi systems to Refs, 
min] in which the phenomenon of the revivals after a quantum quench is construed as a recombination of the 
fastest quasiparticles in the system. Specifically, the speed of the fastest excitations in the system is upper bounded 
by the Lieb-Robinson speed and this upper bound gives a lower bound to the revival time |22j . On the other hand, 

the most distinct revivals are given by the steady values of the group velocity Vg{k) = ED and in 

this scenario the first revival is the one corresponding to the maximum group velocity (vmax = mhXk Vg(k)) and the 
revival time scale could be given by the following estimate 


N 

Trev = -Z -• (9) 

^^max 

To investigate the revival times in LE of the extended quantum compass model in a transverse field, we have plotted 
the LE for various system sizes in Fig. [^for different types of quenches. It is easy to see that increasing the length 
of the chain decreases the oscillation of the LE and the height of first revival decreases gradually and, finally, vanishes 
as Af — 00 . Examining the details shows the scaling behavior of first revival time Trev versus N. This is plotted in 
Fig.[TT] which shows the linear behavior of Trev versus N. The scaling behavior is 
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T’ — KN^ 

rev — -tv 1 V , 


( 10 ) 


with exponent (5=1. 

A similar analysis can be carried out to obtain the K coefficient. Our calculations show that K = — where 

Vmax are given by following estimate 


, diEl±El) 

' dk ' 


.) = 


,dEo 
' dk 


.dEi 

' dk 




In Fig. the first derivative of the ground state Eq and excited state Ei with respect to the momentum k has 
been depicted for different critical points. The numerical calculations show that Vmax = min(| | Imaa:) 

and the revival time shows the universality and scaling around the quantum critical points consistent with those 
predicted by Eq. (10) and Treviho)/Trev{hTr) = VmaxihTr)/Vmax{ho) = 2.61 which consistent with those obtained in 
Fig. 1^ A surprising result occurs in the critical surfaces which are the boundaries between the spin-flop (I) and the 
saturated ferromagnetic (III) phases and between the strip antiferromagnetic (IV) and the saturated ferromagnetic 
(III) phases which correspond to = \/(Ji -I- L 2 ){J 2 — Li). In this critical magnetic field, the revival time of 
the LE is inversely proportional to the the maximum group velocity of the excited state quasiparticles, while in 
ho — yj (Ji + L 2 ){J 2 + Li) critical field the revival time is proportional to the inverse of the maximum group velocity 
of the ground state quasiparticles. It should be pointed out that the regions of phase space with a different Lieb- 
Robinson speed do not correspond to different quantum phases. Quantum phase transition happens in the ground 
state at a singular point of the phase space, whereas the Lieb-Robinson speed is the maximum speed of any signal in 
any state and not just small excitations over the ground state [15] . 


VI. SUMMARY AND CONCLUSIONS 

In this work we have studied the ground state fidelity, fidelity susceptibility, and quench dynamics of the one¬ 
dimensional extended quantum compass model (EQCM) in a transverse field. We use the Jordan-Wigner transforma¬ 
tion to construct an explicit analytic expressions of the ground state fidelity and Loschmidt echo (LE) of this model. 
We show how the fidelity susceptibility and LE could detect the quantum phase transitions in the inhomogeneous 
system. Moreover, we have investigated the universality and scaling properties of fidelity susceptibility and revival 
time. The results show that the fidelity susceptibility exhibits beautiful scaling law close to the critical points with 
exponent u = 1 exactly corresponds to the correlation length exponent of the Ising model in a transverse field. 

We show that the LE exhibits a universal structure of revivals that is independent of the initial state and the size of 
the quench. The information travels through the spin system via wave packets of quasiparticles, and the first revival 
appears when the wave packets propagates with the minimum of the group velocities of ground state and excited 
state. This interpretation explain the universality of the revival structure since group velocities only depends on 
the dispersion relation of the quasiparticles related to the ground state and excited state of the quench Hamiltonian. 
Depending on the critical surfaces, the structure of the revivals after critical quantum quenches represents two different 
equilibration dynamics, whereas examination of the fidelity susceptibility shows the same correlation length exponent 
for all critical points. 


Acknowledgments 

The author would like to thank G. Watanabe, A. Akbari and V. Karimipour for reading the manuscript and valuable 
comments. 


VII. APPENDIX 
A. Unitary Transformation 

The unitary transformation matrix U which can transform the Hamiltonian Eq. ([^into a diagonal form, has the 
following form 
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( 

^‘^,El 

^3, El 

1\ 

Ul,-El 

U 2 -EI 

^3-El 

1 

Ui,El 

^2,El 

^3,El 

1 

V ^1-El 

^2 -El 

U 3 -EI 

1/ 


where 


Ul,±E 


U2,±E^ 


U: 


S,±E- 


2 h±E'^ L* 

J* J* 

L*(2/i±S^) j 

[i2hTE^)\ 

(^L( J*2 - L»2) _ ± £:“)2^ 

)] 


* 

to 

1 

,2 

J*{ja - 

- |L|2) - JL*2 



— L''{ 2 h±E“) 1 1 


- L*"^)-L*{ 2 h±E'^f^ 

)] 


L*- 

j*2 

J*( 7 a- |L| 2 ) - JL »2 



L(J*^-L*2)-L*(2h±S 

f)’); 


J*( 7 a-|L| 2 )-JL *2 



a = q,p, j = 1 for a = q and 7 = — 1 for a = p. 


B. Ground State 

By using the unitary transformation the unnormalized eigenvectors and eigenvalues have the following expressions 
in terms of the coupling constants: 


IV'm) = n 

k 


vm+vT cfc!.\|0) +U 3 ™ C«VJ,|0) +< c«J,c^^|0) +< c^V\|0) +< cfc^^lO) 
vT cf m = 0,1, • • • , 7. 


En 


7)0 

V 4 


-72^4^2 +I J|2 + |L|2 + C, 

64h^ - J*^L^ + J2(2J*2 _ L*^) + 2|JpC - 16Eoh^ - 4h{\J\^ - \L\^ + C)Eo + 8h2(3|Jp - 2|L|2 + 2C') 

8h2| J|2 + 2| jpc + 2| J|4 - J2L*2 _ j* 2^2 
^.o_ Eo{4h-Eo){LJ* - JL*)/2 

8/l2| J|2 + 2| jpc + 2| J|4 - J2L»2 _ j*2^2 ’ 

(4/1 - Eo){4h‘^J + J*( J2 - L2) + JC) 

~8h2| J|2 + 2| jpc + 2| J|4 - J2L*2 - J*22,2 ’ 

(4h - £;o)(4h2 J* + J(J*2 _ L*2^ 

~8h2| J|2 + 2| jpc + 2| J|4 - J2L*2 - J*22,2 ’ 

7)0 =0, = 1, 
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El = -^/2V4/l2 + |J|2 + |L|2_C', 

-64/1^ + J* 27,2 _ j2(2j» 2 _ 2,-2) 2| Jpc + 16Eih^ + 4/i(| Jp - [Lp - C)Ei - 8h^{3\J\^ - 2\L\^ - 2C) 

“ -8/l2| J|2 + 2| jpc - 2| J|4 + J2L*2 j*2^2 

1 1 Ei{Ah- Ei){LJ* - JL*)/2 

^2 - ^5 “ “ -8/i 2| J|2 + 2| jpc - 2| J|4 + J2L*2 j*2^2 ’ 

1 {4h-Eo)i4h^J + J*{J^-L^)-JC) 

~ ~ -8/l2| J|2 + 2| J|2C - 2| J|4 + J 27^»2 j*2^2 ’ 

1 (4/i - £;o)(4/i2 J* + J(J*2 - L*^) - J*^) 

~ -8/l2| J|2 + 2| jpc - 2| J|4 + J22^»2 j*2^2 ’ 

^6 = 1'7 = 0> ^8 = 1. 


E 2 = 0, vl=V2=V2=V2=vl = V2= 0, V 2 = -ff = -1, 

E 3 = 0 , vl=vl = vl = v^ = vl = vl = vl= 0 , vl = 1 , 

E4 = 0, vl = vl = vl = vj = vl = vl = vl = 0, vl = 1, 

E5 = 0, vl = vl = vl=vl = 0, vl = -vl = -1, 

3 _ 4hL 4 _ 4hL* 

~ LJ* + JL* ’ “ LJ* + JL* ’ 


Er 


vl 


y2v/4/i2 + I J|2 + |L|2 - C, 

-64/1^ + J»27,2 _ j 2(2J*2 _ ^*2^ 2| Jpc + + 4/i(| Jp - |L|2 - - Sh^ii\J\^ - 2\L\^ - 2C) 

-8/l2| J|2 + 2| jpc - 2| J|4 + J2L»2 j*2^2 

6^_ £;6(4fe-g6)(4^J*-./4^*)/2 

-8/l2| J|2 + 2| jpc - 2| J|4 + J2L»2 + J*2/^2 ’ 

(4/1 - Ee){4h^J + J*( J2 _ l 2) _ 

-8/l2| J|2 + 2| jpc - 2| J|4 + J2L*2 j*2^2 ’ 

(4/i - Ee){4h^J* + J( J*2 - L*^) - J*C) 

-8/l2| J|2 + 2| jpc - 2| J|4 + J2L*2 j*2^2 ’ 

= 0, t'g = 1, 






y2V4/i" + I^P + |ip + C^ 

64/i4 - J*22,2 + j 2(2 j *2 _ 2,-2) 2|Jpc - l%E^h^ - 4h{\J\^ - \L\^ + C)Ey + 8/72(3|Jp - 2\L\^ + 2C') 

8/l2| J|2 + 2| jpc + 2| J|4 - J2L*2 _ j*2^2 
7 Er{4h- Er){LJ* - JL*)/2 

~ -8/l2| J|2 _ 2| jpc - 2| J|4 + J2L*2 j*2^2 ’ 

(4/7 - Ej){4h‘^J + J*( J2 - L2) + JC) 

~8/72| J|2 + 2| jpc + 2| J|4 - J2L*2 - J*22,2 ’ 

(4/7 - £;7)(4/72 J* + J(J*2 - L*^) + J*^) 

~8/72| J|2 + 2| jpc + 2| J|4 - J2L*2 - ’ 

V7 =0, 7;| = 1, 


in which C = \/a? — b. 
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